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ON ONE DIMENSIONAL INVERSE PROBLEMS ARISING FROM 
POLARIMETRIC MEASUREMENTS OF NEMATIC LIQUID 

CRYSTALS 

YVES CAPDEBOSCQ AND BASANG TSERING-XIAO 



Abstract. We revisit the problem of determining dielectric parameters in 
layered nematic liquid crystals from polarimetric measurements originally in- 
troduced by Lionheart & Newton. After a detailed analysis of the model, 
of the scales involved, and of natural obstacles to the reconstruction of more 
than one dielectric parameters, we produce two simple one-dimensional inverse 
problems which can be studied without any expertise in liquid crystals. We 
then confirm that very little can be recovered about the internal configuration 
of smooth dielectric parameters from these measurements, and give a unique- 
ness result for one of the two problem, when the unknown parameter satisfies 
a monotonicity property. In that case, the available data can be expressed in 
terms of Laplace and Hankel transforms. 
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1. Introduction 

Consider the experiment sketched in Figure [TJ A polarized focused laser beam 
is split in two beams. One beam propagates in an homogeneous medium (in light 
blue), and travels through a slab containing the liquid crystal (LC) cell. The 
transmitted beam exits on the opposite side of the slab, and enters an analyser, 




Figure 1. Sketch of the polarimetric measurement experiment. A 
polarized laser beam is split in two, the main part is focused, sent 
through a thin slab containing the liquid crystal. The slab is placed 
in an isotropic, homogeneous medium whose index is close to the 
average index of the slab. The transmitted beam is refocused, and 
compared to the initial beam by an analyser. The planar cross- 
section sketched above is referred to as the (xi-x^) plane in the 
text. 

l 
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where it is compared to the other part of the split beam, which did not pass through 
the slab. The analyser delivers four real numbers, called the Stokes parameters. 
The LC cell is placed on a cylindrical mount, to allow variations of the incident 
angle of the laser in the slab. The Stokes parameters, appropriately normalized, 
are collected for a large range of incident angles. The purpose of this article is to 
study the dependence of the data on the dielectric tensor of the liquid crystal, and 
discuss its possible reconstruction. 

The liquid crystal contained in the slab has a layered structure: the dielectric 
permittivity is constant and anisotropic in each layer. The magnetic permeability 
is a constant number. The liquid crystal is coated with thin films (polyimidc, 
indium tin oxide...) but these thin hims will not be taken into account in our 
study. This experiment was performed at Hewlett-Packard Laboratories, Bristol 
and was analysed by Lionheart & Newton [Tj using a singular value decomposition 
approach. 

To model this problem we adopt the approximation introduced in Berreman [2] , 
following Lionheart & Newton [Tj . The incident polarized laser beam in the isotropic 
medium in front of the liquid crystal is modelled by monochromatic plane waves 
incident obliquely in the (xi-x 3 ) plane (the plane in which the experiment is 
sketched), with direction s = (sin(0),O,cos(0)). The corresponding electric field 
has two components, 

S\ = £1 exp (iw (— ■? ■ it] 

which is parallel to s]_ — (cos(0), 0, — sin(0)), and 

S\ = E 2 exp (iw (— -f ■ it] 

along s \ = (0,1,0). To take into account possible phase differences, the ampli- 
tudes of E\ and E 2 are complex numbers. When the incident field is polarized lin- 
early, E\ and E2 are of the form Acos((f>) and Asm(<f>), with A £ C and <fi £ [0, 2ir). 
When the incident field has a circular polarization, E2 — ±iE\. The central 
frequency of the laser is denoted by w, c is the speed of light, and no is the index of 
the isotropic medium around the slab on the wave path. We assume that within the 
slab, the propagating electric and magnetic fields £ and T-L stay alike plane waves 
along the x\ direction, and do not depend on xi- Thus £ and T-L are modelled by 
the following ansatz 

(1) £ = E(x 3 ) cxp{iw( — S1X1)}, T-L = H(x 3 ) exp{iw( — sixi)}, 

c c 

and satisfy the time-harmonic Maxwell's equations 

(2) VA£ = icj^ H, V AH = -iuJS Q e(x 3 )£. 

Where e is the relative dielectric tensor, with constant entries outside the slab 
and varying as a function of x 3 within the LC cell, and fio and £0 are universal 
constants. See Lionheart & Newton pQ, Tsering-Xiao [3] and references therein for 
more details and discussions on this model, and Lavrentovich [J for further insight 
on optical measurements of liquid crystals. Substituting ([TJ into ((2J, Maxwell's 
system becomes a 4 x 4 system of ordinary differential equations in x 3 where the 
unknown is the vector X = (f 1 ,H 2 ,^2, ~Tii) T which satisfies 

(3) ^=i-M(x 3 )X(x 3 ). 

ax 3 c 
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The Berreman vector X is the tangential part of the electromagnetic fields and 
represents the optical field. The Berreman matrix M depends on the dielectric 
tensor e. In the outer isotropic medium, the Berreman matrix has a block diagonal 
structure, corresponding to the two possible polarizations of the electric field, along 
i j_ and i \. The Berreman matrix has four eigenvectors (14)i=i,..,4- The vectors 
V\ and V3 correspond to the incoming electric and magnetic fields, whereas the 
vectors V% and V4 correspond to the outgoing electric and magnetic fields, see 
Section 12.11 The direct problem is formulated in Lionheart & Newton [I] as a 
transmission problem. The incident field and the transmitted field are represented 
by the input vector X % = £{V\ + £\V^ and the output vector X 1 — £\V\ + £|^- 
The reflected field X r is a combination of the outgoing eigenvectors V% and V4. If h 
is the thickness of the slab, and P{h) is the transmission matrix coming from (|3|). 
the transmission problem takes the form 

(4) P{h){X l +X r )=X t . 
The Stokes parameters (the measured data) are 

(5) S(E{,Ei9,e) = (\£[\ 2 + |f*| 2 , \£[\ 2 - \S l 2 \ 2 , 2S(£j 

The questions investigated numerically in Lionheart & Newton [I] are 

(1) Can the dielectric tensor through the LC cell be deduced from the data? 

(2) How is the solution affected by the range of incident angles and input 
polarisations used? 

(3) Given the limited accuracy of the polarimeter how much information can 
be deduced from the data? 

The goal of this paper is to address these questions analytically. In section [2] 
we reformulate the problem in an equivalent, non-dimensional form, taking into 
account the scale of the various parameters of the problem. In section |2~T1 we detail 
the derivation of the Berreman model ((3|), given in a non-dimensional form by (|7|). 
We consider two types of LC cells, orthorhombic and nematic uniaxial. In sec- 
tion 12.21 we investigate the case of an orthorhombic medium with principle axes 
aligned with the coordinate axes, that is, we assume that e is a diagonal matrix- 
valued function. This first model was also discussed in Lionheart & Newton [1] , with 
a different approach. Orthorhombic crystals are one of the typical anisotropic ma- 
terials. Within the class of liquid crystals, orthorhombic symmetry was considered 
as a convenient theoretical possibility in early developments, see e.g. Freiser [5]. 
Since then, it has been observed for specific bi-axial liquid crystals by Hegmann et 
al. [7], but it is deemed to be very rare see e.g. Karahaliou et al. [8]. In sec- 
tion 12. 3( we turn to the case of a nematic uniaxial LC cell. This is a very common 
model for LC cells. In that case, e can be expressed in terms of a director profile, 
that is 

(6) e = e±I d + (e|| -£jjn®n, 

where, for each 2:3, 11(2:3) G S 2 is a unit vector, the so-called director vector of 
the LC cell. The two constant eigenvalues of e are usually written e± and en. 
The numbers y / £± and v /eTT are known as the refractive indices of ordinary and 
extraordinary waves. The subscripts || and _L refer to field directions respectively 
parallel and perpendicular to the director vector n. Formula (|6|) implies that the 
electric field energy propagates separately along the ordinary and extraordinary 
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direction. The derivation of this model and further background can be found in 
Virga [9], de Gennes & Prost [10] and Chandrasekhar [TT] . 

For both models, we highlight intrinsic obstacles to the determination of the 
dielectric parameters: non-uniqueness is endemic, regardless of practical limitations 
such as the range and precision of the measurements involved. In section I2.4[ we 
express the scales of the various coefficients involved in the experiment described 
in Lionheart & Newton [1] in terms of a non-dimensional parameter. 

We conclude this section by producing three one dimensional inverse problems 
coming from both the orthorhombic LC cell and nematic uniaxial LC cell model, 
Problem A given by (|29[) . Problem B given by (|30l) and Problem C given by (f3Tj) . 
Problem A and Problem B concern the reconstruction of one coefficient from the 
amplitude of the transmission data. This coefficient would be £22 for the orthorhom- 
bic LC cell example, and one component of director profile n a priori assumed to 
lie in the (xi-xs) incident plane. Problem C presents the independent phase data 
available from the nematic uniaxial LC cell model, which is just one numerical 
value. 

Section [3] is devoted to Problem A. Under a smoothness assumption of the co- 
efficient to be reconstructed, we describe what can be extracted from the data 
corresponding to moderate angles of incidence. We then explain how Problem B 
connects to Problem A, as a first order approximation, when the full range of angles 
of incidence are considered. 

Section S] is devoted to Problem B. For this model problem, we can describe the 
class of equivalent parameters, and give a uniqueness result with a monotonicity 
assumption. We show that this problem can be reformulated in terms of the classical 
Laplace and Hankel transforms. 

Section [5] contains the proofs of various technical intermediate results given in 
the previous sections. We summarize our findings and discuss possible extensions in 
section [6] We use some facts concerning systems of ordinary differential equations 
in this paper, they are given in appendix for the reader's convenience. 

2. Modelling and scaling assumptions 

2.1. An alternative equivalent formulation. In this section, we prove the fol- 
lowing proposition. 

Proposition 1. The transmission problem ([3" T) -(|4 | can be written in a non-dimensional 
form as follows. Given I\ and I 2 in C, find T\ 1 T 2l R\ and R 2 in C such that there 
exists a continuous solution of 

1 ' ry -1 

— = i-B(t)Z(t), te[o,i], 

at i] 

(7) Z(0) = [I 1 -R l ,I 1 + R l ,I 2 -R 2 ,I 2 + R 2 ] T , 

Z{1) = [T u T u T 2 ,T 2 f , 

where B is a matrix-valued piecewise continuous function of t € [0, 1] given by 
(8) 

-sin(0)£ 13 A " + £ ~ 33 = 1 - sin^A" 1 ^ 



1 
B = — x 
£33 



A (en£33 - £13) - sm(6>)ei 3 £12233 - £13223 

S S S S m-n(ff\\-^-S H (A + £~22 ~ l)£~33 ~ £23 

£l2£33 — £l3£23 —Sm(0)A £23 U 3; 

A£ 33 
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with 
(9) 



h' 



1 



e = 



Oo) 



r£> 



cun h 



= 7] , and cos(#) = A. 



T7ie length h is the thickness of the LC slab, 9 is the angle of incidence, and n is 
the index of the homogeneous isotropic medium surrounding the slab. 

The Stokes parameter data given by ([5]) is equivalent to 

(10) S{I u I 2 ,e) = (7\ exp (i/(Ii, I 2j 0)) , T 2 exp (i/(Ii,I 2 , 0))), 
where the real function f is unknown, and may depend on I\,I<i and 9. 

Problem ([7]) has only partial conditions on the solution at the start point and 
partial conditions at the end point. It is therefore not a standard initial eigenvalue 
problem for system of ordinary differential equations (|14[) . It is nevertheless well- 
posed, as the following proposition shows, proved in \K^ for completeness. 

Proposition 2. Given I\ and I<i, there exist a unique pair of reflection parameters 
Ri and R2, a unique pair of transmission parameters T\ and T%, and a unique Z 
solution of (f7|). 

Proof of proposition^ We rescale the space variable x% by t = x$/h. Problem [3] 
becomes 

(11) ^ = -M(t)X(t), 

at r\ 



and M is given by 
(12) 



M 



"o 



»o 



E13 sin(g) 

£33 
eoc(en£33 — e? 3 ) 

£33 



^QC (£12 £33— £13 £23) 



£33— sin(fl) 

n £0C£33 

£13 sin(fl) 
£33 



£23 sin(fl) 



£23 sin((?) 



£33 

£QC(£12£33~£13£23) 

£33 





n 



"0 

£0 C (£22£33 



-sin(e) 2 e 3 3) 



In the outer isotropic medium, the matrix M simplifies to 



Mi 





n e c 







«0£0C 





-A 









uqEocX 



n soc 





which has a block diagonal structure, corresponding to the independence of two 
possible polarisations of the electric field. The matrix Mi has four eigenvectors, two 
corresponding to the incoming electric and magnetic fields, and two corresponding 
to the outgoing electric and magnetic fields. Ordering them in agreement with the 
block structure, these are 



V 



corresponding to the eigenvalues A, —A, A, —A. We note that 
X 1 = £{ V [1, 0, 0, 0] T + El V [0, 0, 1, 0] T . 



A 


-A 








e cn 




e cn 






1 

£qCUq\ 




1 

— eocngA 
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The transmission problem takes a simpler form if the matrix M is written in the 
eigenbasis of M\, that is, 

(13) M = V _1 MV. 
The transmission problem ©-(Hi then becomes 

Given an initial condition I = [I\, 0, I2, 0] , find T — [Ti, 0, T2, 0] 
and R — [0, i?i,0, -R2] suc/i t/iat i/ie solution of 

(14) ^ = iM (t)Y(t)) te[0 ,i], 

satisfies 

Y{0) = I + R, Y(1)=T. 

The stokes parameter are then given by 

s{h,i 2 ,e)= (|T 1 | 2 + |r 2 | 2 ,|r 1 | 2 -|r 2 | 2 ,2T 1 ^ 

A computation shows that the available data is equivalent to the knowledge of the 
vector S given by (|10|) . Finally, the structure of the matrix M, is simpler after a 
block — 7r/4 a rotation. We define 



where 



B = 


R^(M)R_ i 




1 

V2 


"1-10 
1 1 
1 
1 






-1 
1 



R2L = 



and the entries of B are given by (|5|). After these simplifications, problem (fT4")l 
transforms into problem ([?])• □ 

2.2. Orthorhombic LC cell model. The following proposition summarizes our 
findings concerning the orthorhombic model. 

Proposition 3. If e is a diagonal matrix with three independent piecewise contin- 
uous functions (£u)i=i,2,3 as diagonal entries, not all three can be determined from 
the Stokes parameters: en and £33 cannot be determined independently from each 
other. On the other hand, the transmission coefficient T\ is independent of £22, 
whereas T 2 is uniquely determined by the system 

V X^- = i(A 2 + £ 22 - l)z 4 in [0,1], 

(15) r,^ - iA*3«i[0,l], 

at 



z 3 (0) - I-R2, z 4 (0) = l + i? 2 , z 3 (l) = z 4 (l)=T< 



2- 



If en and £33 are unknown, the available data to determine £22 is IT2I; where \ ■ \ 
denotes the complex modulus. If the phase ofT\ is known, and this is the case when 
both £n and £33 are known, the available data is T2. 
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Proof. If we assume that e is a diagonal matrix-valued function with en > C > 
for i = 1, 2, 3 and C is a positive constant, the Berreman matrix B becomes block 
diagonal: 

A 2 + £33 - 1 



B = 




Affn 






£33 A 











£22 



By inspection we note that the unknowns (i?i,Ti) only depend on I\, and the 
unknowns (i?2,72) only depend on I<i'- the two transmission modes are decoupled. 
By linearity, we can thus set I\ = 1% = \. The transmission problem relating R 2 
and T 2 becomes (fT5|) . 

To write the transmission problem relating R\ and T\ in a similar form, we write 
Ai = J Q Eii(s)ds 7 and introduce the change of variable Ai^(x) = J Q Su(s)ds, and 
write the inverse change of variable [i = v~ x . The first block of ([7]) then becomes 



(16) ^ 



■nX dz\ 
Ai dv 

r\ dz 2 
Ai dv 

21 (0) 



. A 2 +£33(/i(^))-l 

eu(MM)e33(^M) 
iAzi in [0,1], 
1--Rl, Jf2(0) = l 



z 2 m 



[0,1], 



R x , z 1 (l) = z 2 (l) = T l . 



It is clear that problem (|16|) contains too many unknown functions and parameters 
for each of them to be uniquely determined by the map A —¥ T\. Furthermore, 
formula (flu]) given in proposition Q] shows that for each A the pair (Ti,T2) is only 
known up to an arbitrary phase shift. Since T\ and T2 depend on different unknown 
functions, in general only |Ti| and IT2I are available. □ 

Proposition 4. The transmission coefficient Ti is unchanged if t —} £22 (t) is 
replaced by t —t £22(1 — t). 

We check this classical reversibility property in [B] 



2.3. UniELxial nematic LC cell model. Let us now consider the nematic uniaxial 
model ([5]) . The index of the surrounding homogeneous medium is given by 



(17) n 2 a = 

We parametrize the director vector by 



V^l£|| 



n = (cos(^) cos(y), cos(V') sin(^), sin^)), 



(18) 



(V>,v)e (- 



2'2 J 



[0,77). 



The tilt angle ifi and the azimuthal angle ip parametrize a half-sphere only, since 
e is unchanged if n is changed into — n. Uniaxial configurations are a priori less 
complex than the orthorhombic ones, since two functions instead of three are to be 
determined. We simplify the problem even further, and impose that ip — 0: the 
director vector n lies in the incident plane. This is still not sufficient for uniqueness, 
as shown by the following proposition. 
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Proposition 5. If e is the dielectric permittivity matrix of a nematic uniaxial 
LC cell modelled by ((6|), the Stokes parameter data are not sufficient to determine 
uniquely the tilt angle ip and the azimuthal angle if of the director vector defined 
in (|18|) . When the director vector stays in the incident plane, that is <p = 0, the 
available data are 



where 
(19) 



T(A)F(A) and T(X)/F(\), 

sm(2ip(s)) 



F(X) = exp M?rWl - A 2 I — 



-ds 



<5cos(2V>(s)) 
The parameters with A, and r\ are defined in ©, whereas S is given by 

£j_ — £|| 



e±+£|| 




The map A — > T(A) is defined by 




r^ i(V 1 S (nr(««lt)) 1 S 




' ' dt l \y ' vi-^ 2 V [ { '' i + VT- 


-5 2 


duo 
(21) fj— - = iXui in [0,1], 





U2 m 



[0,1], 



ui(0) = l--R(A), u 2 (0) = l + R(X), ui(l) = «2(l) = T(A). 

where /«(t(£)) = £ /or a^ £ G [0, 1], u/ii/i 



r(i) = 



1 



-ds, A/" = 



-ds, 77 = 



n 



Af J I + 5 cos(2ip(s)) " ' J 1 + 6 coa(2ip(s)Y 

The map t —¥ ip(t) is not uniquely determined by the available data. 
Proof. When ip = 0, the Berreman matrix B is block diagonal, with 

s/T~P{\ 2 -l)\ 



VT^PAT 



(B)i 



and 



(B) 



5sin(0)sin(2V>) 1 / 

<5cos(2^) + l A^ 1 ^ o~cos(2i/>) + l 

\/l-6 2 \ * S sin(6») sin(2^) 

<5cos(2'0) + l (5cos(2-0) + l 



3. .4,3. .4 





A 



e_L 



n 2 A 





The second transmission parameter, T 2 , can be computed explicitly independently 
of the tilt angle ip. Thus, using formula (flQ|) we see that the available data to 
determine ip is Xi, and not just its modulus. Performing the change of unknown 



(ui{t),u 2 {t)) = exp 



,sin(6») 



<Jsin(2^(s)) 



77 Jo o~cos(2i/>(s)) + 1 



ds (z 1 (t),z 2 (t)), 
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Figure 2. Example of non- uniqueness for nematic crystals. If 
the azimuthal angle ip is constantly zero, then Ti(ip\) — Ti{_ip2) for 
any incident angle 9 . 



we see that I\ is given by T(X)F(X) or T(X)/F(X) with T(X) determined by 



^ = ■ ^ 



(22) 



V 



du 2 



iX- 



y/l-S 2 



l + (5cos(2V>) 
-Ui in [0, 1], 



u 2 in [0,1], 



dt *"l + <5cos(2^) 

mi(0) = l--R(A), u 2 (0) = l + i?(A), u 1 (l) = « a (l) = r(A). 



Since both positive and negative incident angles are measured, but T(A) only de- 
pends on A = cos(0), the data T x is both T{X)F{X) and T(A)/F(A). 

We then change t to r in (|22|) and obtain (|2Tj) (where we named t instead r the 
dummy variable). 

The coefficient T(A) depends on ip, through cos(2^/;) : positive and negative tilt 
angles are not distinguishable. The phase difference between T(A) and T(X)F(X) 
or T(X)/F(X) depends on -0 via the number 



K = 



sin(2tp(s)) 



-ds. 



/o l + 5cos(2^(s)) 

In Figure[2]we show the graph of two tilt angle functions which would have identical 
transmission parameter T(X) for all A. To construct these examples, we chose to 
repeat a given compactly supported pattern five times. In two instances out of five, 
the pattern is flipped with respect to the t axis. There are ten possibilities: we 
chose two different ones arbitrarily. The parameter T(A) is equal for both graphs, 
as it is invariant under arbitrary changes of the tilt angle sign. The constant K 
sums the contribution of the five patterns, three with a plus sign and two with a 
minus sign. The way these patterns are ordered does not change this integral. This 
non-uniqueness comes in addition to the one already highlighted in proposition 31 
namely that T(A) is unchanged if t — > ip(t) is replaced by t — > ip(l — t), see|B] □ 



2.4. Scaling assumptions. In this section, we discuss the scales of the various 
quantities involved. The numerical values of the physical parameters used in this 
section are taken from Lionheart & Newton pQ. 

Frequency. For a He-Ne laser of wavelength 0.633 /jm and a slab of thickness 
5 /xm, we find 77 = rj c w 8.4 x 10~ 2 . In what follows, we will therefore assume that 



(23) 



77 is a small dimensionless parameter. 
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Dielectric parameters. In the case of a orthorhombic (diagonal) dielectric 
tensor, we will assume that 

(24) e 2 2 = 1 + T] a q, where sup \q(t)\ < 1, and a near 1. 

te[o,i] 

A similar assumption for the nematic uniaxial LC cell to give a simple form to ()19j) 
is 

(25) 5 = r] a . 

For problem (12T1) to match the orthorhombic assumption (|24|) , we can choose 

For the nematic uniaxial LC cell the values are | ^/ejT — a/^TTI ~ 0-15, n o ~ 1.52. The 
scaling ([25]) yields a ps 0.94 whereas (|26|) leads to a ~ 0.92, both values are indeed 
close to one. 

Incident angle. We write the range of the incident angle as follows: 

(27) A G [■s/rrf', 1], and r > 1. 

In the experiment considered, varies between from 0° to about 70° : for larger 
angles the measurements become unreliable. The extremal value (cos(70°)) 2 corre- 
sponds to r = t c s=a 1.2. 

Measurement error. We will assume that the measured data is accurate up 
to errors of order 

(28) O (^ 

Where for any x > 0, |C(a;)| < Cx, where C is a constant independent of a;. In the 
experimental case considered, this corresponds to a precision of the order of 10~ 7 
for a normal incidence, and of the order of 10 -4 for the most slanted incidence. 
This assumption models the fact that the measurements become less accurate as 
slant of the slab increases. 

Sampling rate. We suppose that A is measured with a fine sampling rate, 
e.g. T] 2 . Experimentally, 200 incident angles 6 between —70° and 70° are used, 
corresponding to sampling rate 77 1 ' 86 . 

2.5. Model Problems. To summarize the discussion of the previous section, we 
now write down three traceable reconstruction problems pertaining to the or- 
thorhombic LC cell model and the nematic uniaxial LC cell model. 

Problem A. Let r > 1 be a constant, a a parameter close to 1, 
and ?)>0a small parameter. Let q : [0, 1] — >• [— 1, 1] be a piecewise 
continuous function. For every A G [y/Tif, 1], let T(A) and R(X) 
be the solutions of 

V \^- = i(A 2 + ?7 ^) U2 in[0,l] 

(29) 77^ - iAuiin[0,l] 

at 

ui(O) = I-R2, u 2 (0) = l + R 2 , u 1 (l)=« 1 (l)=T(A). 
What can be determined about q from V A { A) = |T(A)|+0 (?7 5 + Q A- 5 )? 
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Remark 6. We note that the lower bound on A, coming from experimental con- 
siderations, has a natural analytic interpretation. When A 2 3> i] a q, Problem (j29|) 
corresponds to a wave propagation problem, slightly perturbed by q. We are there- 
fore in an 'optical' regime: the transmitted wave is very similar to the incident 
wave. When A 2 and Vj a q are of the same order: the parameter q is no longer a 
perturbation but a leading order term. Problem (|29j) becomes diffusive, and one 
can therefore expect the polarimetric measurements to become unreliable. 

Problem A is directly inspired by proposition [3] for the orthorhombic LC cell 
model. We argued that when the dielectric tensor is diagonal, with principle axes 
aligned with the coordinate axes, we cannot hope to reconstruct all three diago- 
nal coefficients. In particular, the entries E\\ and £33 cannot be determined inde- 
pendently. According to proposition [3] the second diagonal entry £22 determines 
uniquely |T 2 |, independently of e\\ and e 33 . The scaling assumption were discussed 
in section 12.41 

Problem A is also relevant for the in-plane nematic uniaxial LC cell model, 
that is when n = (cos(V'), O,sin(-0)). In that case q — cos(2t/>). More precisely, 
proposition \5\ shows that |T(A)| = \T(\)F(X)\ is measurable, and determined by 
Problem A using the scaling assumption discussed in section l2~4l with fj instead of 
i] and 

q = ( cos(2V>(m)) + : |=s? ] (1 



1 + s/T~P J \ 1 + sTT^P, 

with S and fj related by (|26l) . and /1 given by proposition [5j At first order, /i is the 
identity and 8 ~ -q a . Problem (|26p is therefore a close variant of Problem (|2^|) . 

We already know from of the various non-uniqueness examples presented in sec- 
tion 12.21 and 12.31 that little can be determined about the interior values of q from 
Problem A. We discuss it further, assuming q is smooth in section [3J We will see 
in section [4] that we can give a more precise answer to the related problem 
Problem B. Let r > 1, a > be constants, 77 > a small dimen- 
sionless parameter, and q a piecewise continuous function such that 
I q I < 1. For every A g [y/rrj^, 1], let T>B[q] be given by 

(30) A -+ J \jl + ^q{s)d S . 

What are sufficient conditions so that T>B[qi] — X'sfe] implies 
qi = <?2? 
Finally, let us consider the phase information available from the nematic uniaxial 
LC cell data. 

Problem C. Let T>c[tp] be given by 

(31) Vcty] - f 1 .T^m * + 0{5% 

J l + 6costp(s) 

where S is a small parameter. Assuming that tp : [0,1] — > [— 7r,7r] 

changes sign once, at s* 6 (0, 1) and \ip\ is known, find possible 

values for s*. 

As both T(X)F(\) and T(X)/F(X) are measured, F 2 (A) is available. This data is 

measured precisely for close to normal incident angles, that is, ■q~ 1 5y / l — A 2 small, 
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we can extract from (|19| the constant depending on ip, which is T>c['2tp]- This 
problem is naturally very simple to solve: we highlight it here to point out how this 
information can be extracted from the data. On the other hand, not much more 
can be obtained from this problem, since the data in this case is just one value. 
Define 

G{t)= r ywi ds . 

Since < IVKs)! < 7r, G is non decreasing, and strictly increasing if ip does not 
equal or n on a set of positive measure. There are two possibilities, depending on 
whether -0 is positive and then negative, or negative and then positive. 

G(s ) = l -G{\) + \{Vc\1>] + <D(S 2 )) or G( Sl ) = 1 -G{1) - \v c \^\ + 0(5 2 ). 

If 6 is small enough so that —(5(1) < T>c[ip] + 0(5 2 ) < G(l) both si and sq always 
exists, are unique if G is strictly increasing. 

3. On Problem A under a smoothness assumption. 

mln this section, we show that if q smooth, namely C 4 ([0, 1]), moderate angles 
of incidence provide information about the endpoint values of q only. The internal 
configuration of q is not decidable as it is stable under suitable re-arrangements. 
We then explain how Problem A leads to Problem B. 

Our strategy is to find an explicit approximate formula for 

(32) -D'a = ^p - 2- 

Notation. In this section, given x > and y > 0, 0(x y ) means \(D(x v )x~ v \ < K, 
where if is a constant depending on ||<7||c 4 ([o,i]): T given by ([27)) and y only. 

Proposition 7. Assume that q G C 4 ([0,1]). There exist r/ > such that for 
all < i] < r/ and < a < 2, Problem A given by (|29p is equivalent to the 
reconstruction of q from 

p» »i-(5^a))" + (5tw)" + ta(i»- + (*(i))- +0 (^ 

where v\ and V2 are given by 



*<" " ^^l CMd8 ) + T^^cjb>"^ 



nm = V c(A,0) 1 C(A.<) Sil 'Ql C(A '' ) '" 



with 

(34) A=i [ l + l^q) , 5=--A 3 — , andC=(A(l + ^B)) ? 



-V4 i fli 2 

A ,, y , B = --A* — , andC=(A(l + ^ 

m > max{j) : 1 + ?7 2 " Q S > for all t G [0, 1]}. 



Remark 8. This is an approximation of WKB-type. The uniform pointwise error 
estimate involves second order derivatives of V\ and v 2 and requires g€ C 4 ([0, 1]). 
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Figure 3. Example of non- unique determination from the data. 
If the dielectric permittivity q follows the blue curve, or the red 
variant, the available data T>' A is the same, for any A > 0. 



We will prove this proposition in section [5] Note that because up to error 
terms, V' A depends only on q via the value of (C, ^) at t = and t = 1 and on 

L C(A, s)ds, only a large class of equivalent q can be determined in Problem A, at 
best. 

Corollary 9. The usable data contains only some information of the first three 
derivatives of q at t = and t = 1. and on 



(35) 



F 



r] a r\ a dq r] a d q 
A 2 " 9 ' A 2 ^' A 2 dt 2 ' 



dt. 



where F is an explicit algebraic function, independent of A. Since (|35[) is stable 
under sufficiently smooth rearrangements of q (see e.g. Example \lCfy the determi- 
nation of q from the data is not possible without additional a priori information on 
the variations of q. 

Example 10. Let qo £ C ([0, 1]) be a positive function which is constant on the 
interval [to,ti] with 1 < to < i\ < 0. Let q\ £ C 4 (K) be positive function with 
compact support in [0, (t\ — £q)/2]. Consider the family of functions q s given by 

q s (t) = q Q (t) for t £ [0,t Q ] U [ti, 1], and q s {t) = q Q {t) + q x {t - s) for t £ (t ,<l) 

where s is a parameter in (to, (to + ti)/2). Then for all such s, q s £ C 4 ([0, 1]) has 
the same end-point values, and for any smoooth functional F, 



ds 



-r F \ T¥^>T¥-7T> 



rf dq s rj a d 2 q s 



A 2 dt ' A 2 dt 2 



dt = 0. 



Figure [3] represents such a construction for two different values of s. Note that the 
order of accuracy is not the issue when q is smooth, as the WKB-type ansatz can 
be continued to obtain any order of accuracy, for a suitable 770- Counter examples 
to uniqueness such as the one depicted in Figure [3] would still hold. 
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The main result in this section is devoted to the extraction of a few properties 
of q at t = and t = 1 from moderate angles of incidence, for a sufficiently smooth 

q- 

Proposition 11. Suppose that q £ C 4 ([0,1]) 7 and that a is close to one. Using 
the data available for \ < A < 1, in Problem A given by (]2T)1) . the parameters 
Ai,i = 1, . . . , 5 given by 

A x = q(0)+q(l) + O( V 2 ), 
A 2 = q(0)q(l) +0( V 2 ), 



A4, = f q{s)ds + 0(r ] 3 
Jo 



A 5 = / q 2 (s)ds + 0(n 4 - 2a ). 
Jo 

can be determined. 

The proof of this proposition is given in section [5j 

Remark 12. This hrst order data respects the invariance by mirror symmetry dis- 
cussed in Proposition |U If the LC cell is flipped, its index becomes q(t) = q(l —t), 
and the parameters A;, i = 1, . . . , 5 are unchanged. The fact that only a handful 
of moments can be recovered from moderate incidence is known for related more 
general problems, see Sharafutdinov [12] . 

To investigate the problem further, let us assume that the map A — > C(X,t) 
defined in ([34]) is known at the endpoints t — and t = 1 up to the order of 
approximation. The usable data T>' A takes the form 

'2A f 1 „ \ .. „ {rf+^ 



V' A = Ki cos ( — J C{\, s)ds + K 2 J + K 3 + O 



A 5 

where K\, K 2 and K3 are known parameters. Since the data is sampled in A at a 
fine T] 2 rate, the leading order A dependent data is 

(36) f C{X,s)ds 

Jo 

which, at leading order, is 2?s[g](A) given by (|3"0)) . Note that this problem also 

arises naturally if T, and not simply its amplitude, was considered in Problem A, 

as (f31)f is precisely the rate of change of the phase of T. 

4. On Problem B : uniqueness under a monotonicity assumption 

We now turn to the reconstruction of q from T>b [q] ■ 

Theorem 13. Consider A — ¥ X>s[g](A) for all A e [y/rfj^, 1] as defined in Problem B 
given by ()30|) . 

• The map q —¥ T>g [q] is invariant under level-set preserving rearrangements 

of q. That is, given two functions qi and q 2 defined on [0,1] such that 

sup|gj| < 1 which satisfy 

meas{qi{t) =y,te [0,1]} = meas{q 2 (t) =y,t€ [0,1]} My £ [-1,1], 
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we have T> B [q\] = T> B [q 2 ]. 
• Ifq G C 1 ([0,1]) is strictly increasing, sup|g| < 1, andq^ 1 € C 1 ([q(0) , q(l))) , 
then T>b [q] determines q uniquely. The reconstruction problem can be for- 
mulated in terms of classical transforms in this case. There holds 

(37) \v B { q }(X) + ^(X) = 4=CF oQ ( X 



A " ,y ' d\ w ^f u ^ V07 5 

where L is the Laplace transform, Cf(s) — L e~ st f{t)dt, and Fq is the 
first Hankel transform, Fof(s) — L tf(t)Jo(st)dt, where Jq is the Bessel 
function of the first kind of order 0, and Q is given by 



•v-i 



Q := x -> l {q{0)Ml)) (x 2 )-j—(x 2 ). 

Remark 14. The Hankel Transform is self-invertible, and stable. Inverting the 
Laplace transform is well known to be an exponentially ill-posed task see e.g. Ep- 
stein & Schotland [13]. Even more so when, as it is the case here, the Laplace data 
is limited to an interval {\/r,ri^ a / 2 ) located away from zero. The data 2?e[g] can 
be written as a linear integral operator acting on -| — in the form 



. dq 1 ,rj a , 



9(1) 



V B (\) = K-±— (-L-) with Kf(x) = / y/T+xyf(y)dy for 0<x<l. 

ay ^ 7 9 (o) 

It is clear that K : C°([q(0), <?(1)]) — > C°([?7 q ,t~ 1 ]) is a compact operator : its 
inverse is therefore unbounded. Including grazing angles measures, corresponding 
to r = 1 would still lead to an unstable problem. 

To prove this result, we note that the available data takes a much simpler form 
after an integral transformation. 

Lemma 15. Given a function q such that defined on [0, 1] such that s\xp xG i -n \q(t)\ < 

1, we can define T> B [q\ as the complex valued holomorphic extension ofT> B [q] on 
the annulus \Jr) a T < |£| < 1. 

Let K, represent the holomorphic extension of 

- f^ 2 

2y^e x / e -t dt 
Jo 



1C:x-+l 



to the same annulus. 

Let V' B [q] G C°°(R) be given by 

Z>b [?](*) = / e~ itq ( s Us for all t G K. 

7o 

Then, for any tg R and p £ (V r7 ? Q > 1); we have 

V B [q](t) = ± J * T> B [q](pe is ) K (ir,- a pHe i2s ) ds. 

The correspondence between T> B [q] and T>' B [q] is one to one. 

We prove this result in section [5] This in turn shows that T> B [q] determines q 
uniquely within the class of smooth monotone functions, albeit in a very unstable 
manner. 
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Proof of theorem \13\ . The first statement is a straightforward consequence of the 
coarea formula applied to the map t — > T>' B (i). The second statement follows from 
the change of variable formula as follows 

^b [«](*)= / e~ H ^d S = f e - H yi iq{QU(1)) (y)^(y)dy, 

Jo Jm a V 

Thus 2?b[?] i s the Fourier transform of -| — on the range of q. Applying the 
inverse Fourier transform we therefore recover the values of q(0), q(l) and -| — , 
and therefore q. We do not claim that the requirement q^ 1 6 C 1 ([q(0) , q(l)}) is 
sharp. 

Let us now turn to the equation satisfied by T> B [q\. We write, 

-V B [q A +— fr 1 A = / ^=rfy. 

A d\ J Q y/lf // " \2 ~ 

Using the same change of variable as above, it follows that 



a " mv ' dx w vr 7o /7TT 1 



2 



T ./o 

2 CF f ' A 



e \7^ F f(r)dr 



as announced. □ 

5. Proofs for Proposition [71 Proposition [TT1 and Lemma [T5l 

We start by expressing the available data in terms of the fundamental solutions 
of the ordinary differential system (j2"5)) . Let ui and u 2 be the two solution of 

d u ' 
(38) i 1 2 ^ + (X 2 + ifq)u l =0, t = l,2, 

, , diii , „ , . duo , , 

ui(0) = l, -^(0)=0, ua(0)=0, »?-^(0) = A. 

We have the following identity 

Proposition 16. . The transmission parameter T 2 determined by (129p satisfies 

2 . 77 dui , , ?7 (fit? / x 

where U\ and u 2 are determined by (|38p . Furthermore, 
4 / 77 du\ . . \ ( rj du 



This is a simple calculation performed in |B| 

Proof of Proposition H] We introduce a perturbed problem closely related to the 
original problem when g is smooth. We note C given by (|34p is well defined and 

2 

bounded provided rj is small enough, namely when 1 + Q*B is bounded below by 
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a positive constant. Assume that r\ < 770, with 770 chosen so that C is bounded. A 
tedious but straightforward computation shows that V\ and v 2 are solutions of 

(39) r, 2 ^ + (X 2 + r, a q)v i = 0, i = l,2, 

«i(0) = l, ^(0)=0, « 2 (0) = > 77^(0) - A, 
with 

(40) ll«-?lloo([o,i])<^, 

where K depends on ||g||c 4 ([o.i])- Introducing the function G(x,y) € C 4 ([0,1] x 
[0,1]) given by 

G(x,y) = (v2{x)v!{y) - v 1 (x)v 2 {yj) {q{y) - q(yj) , 

= ^C(X,x)C(X,y) S[n {~ V I C(A ' S)ds ) ^ - q{y)) > 

we hnd, comparing f|38|) and (|39|) and using Duhamel's Formula, that for i = 1,2 
and is [0,1] there holds 

u*(t) - Ui(t) — / G(i, y) («i(y) - v,(y)) dy = —— / G(t, y)vi(y)dy. 

And from the bound (|40|) we deduce that 

A ~~dt ~ ~dt 



+ \\Ui-Vi\\co(lO,l]) < K- 



c°([o,i]) ^ 5 



Proposition [16] then shows that 



2 . 77 dTji , , 77 dv2 , , , . , , „ / 77 +Q 

* = H-z-0) + T-jrW + ^(!) - Ml) + o ' ' 



T A dt w A dt w w w V A5 

and since T ss 1, this approximation also gives an approximation for T of the same 
order, which is the precision up to which T is known. This in turn leads to the 
formula for T>' A announced in Proposition [7] □ 

To prove proposition [TTJ we focus on what can be recovered from T>' A around the 
normal incidence. 

Proposition 17. Suppose that ^ < A < 1, and a is close to one. Then, the 
parameters Ai,i = 1, . . . , 5 given by 

A x = q(0) + q(l) + O( V 2 ), 
A 2 = q(0)q(l) +0( V 2 ), 

d In a , ., . d In a , . _ . . 

a 3 = Y (1) ~^r (0) + o(7?) ' 

A 4 = f q(s)ds + 0(r 1 3 ' a ), 
Jo 

A 5 = [ q 2 (s)ds + 0(n 4 ~ 2a ). 
Jo 

can be extracted from the data T>' A in this moderate range of incidences. 
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Proof. A Taylor expansion of V' A given (|33|) around r\ = 0, using the fact that a > 1 
shows that 

,2o 



2^(A) = 2+|^^(l) 2 + g(0) 2 -(g(l) 3 +g(0) 3 )^ 



+ ?(0)?(1)^ ((g(0) + ? (1))£ - 2^ cos(2^ + 0(A)) 

|(l) ff (0) - §(0)g(l)) ^ sin(2^ + 0(A)) + 0(^+ 2 "), 

where 

A Z" 1 
0(A) =2-(/ C(A,s)ds-l). 

»7 Jo 
Using the formula for C(A, s) given by (f34f we find that for A e [1 — 7r7V?7, 1], for a 
given N , we have 

0(A) - 0(1) = ^V" 1 f q(s)ds + ^-V 2a - 2 j q 2 (s)ds + G(ry 2 ). 

Without changing the order of approximation, since 77 and a are known parameters, 
the available data thus becomes 



A 2 - 2A 2 - Ax (Af - 3A 



■2 



>f 



'7° ol„„lo^ , 1 ~~ *_<*-! /i , ^ _ ^ 2<> 2 



A 2 \Ai-fr 2) cos ^2- + —rT^ + -arn^'A* 

A 2 A 3 1L sin ( 2 * + ^±t,«- 1 JU + *^±r, 2 <*- 2 A 5 ) + Otf), 



A V V A ' * 4A 3 

The conclusion follows, as we have 0(ij~ 2 ) data points available for T>2- We could 
recover these five parameters either by a least square fit, or by deriving explicit 
formulas. □ 

Finally, we prove lemma 1151 

Proof of lemma \15[ Since \rj a q(t)/X 2 \ < r < 1 for every t £ [0,1], we can write 
using the binomial formula 

^ V^" ! A io 

Introducing the holomorphic function V^q of the complex variable £ on the annulus 

\J"q a T < ICI < 1 given by 

we observe that this quantity is well defined, as the series is absolutely convergent 
for \C\ > y/rj™- As A varies in [^/Tr] a , 1], 2?s[g](A) determines Vb[q\ on a non-empty 
real interval and therefore also on the full annulus. We may therefore define ©s[g] 
as the holomorphic extension of Z>s[g](A). 
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We compute, using the residue formula on the ring of radius ^Jrr\ a < p < 1, that 
for any (St, 



/ (qisTds = — \ P > / V B [q] (pe ls ) e 2ms ds. 
Jo 2^/irT (n+ =) J 



2 

This in turn shows, using the fact that all sums are well behaved, that 

(-k] a t) n f 1 



[ exp(-ir) a tq( S ))ds = V Lil^ / ( Q ( s ))» ds 



±- J " V B [q](pJ s )iC{\p 2 te 2is )d Sl 



where 



£(0 = v^£ 



„=o r (» + i) 

It turns out that /C has a closed form expression on the positive real line, in terms 
of the Error Function. For any x > 0, 

er l dt. 



K(x) = l + ^TTxe x eriy/x = 1 + 2^Jxe x I 

Jo 



The converse transformation, using the same steps backwards, gives T>s[q\ in terms 
oiV B [q]. D 

6. Conclusion 

We have analyzed the polarimetric measurement with variable incident angle 
experiment described in Lionheart & Newton [T] . After a detailed inspection of the 
various scales involved in the problem, we produced two related one-dimensional 
reconstruction problems, Problem A given by (|29|) and Problem B given by ([30]) . 
We argued that these problems are relevant for the two Liquid Crystal Cell config- 
urations we considered, orthorhombic and nematic uniaxial. Both problems have a 
simple formulation, and can be studied by applied analysts with no prior exposure 
to Liquid Crystals. 

We partially addressed Problem A, under a smoothness assumption. Further 
investigation of this problem would be of practical and theoretical interest. In the 
smooth case, it would be interesting to see how the boundary data can be recovered 
in an optimal manner numerically. The non-smooth case is left open. While we do 
not believe that much more can be done with very general parameters, it is very 
possible the ill-posedness of this problem is dramatically reduced when additional 
constraints are imposed, e.g. when the medium is piecewise linear. 

We addressed Problem B, by characterising a class of equivalent parameters, 
and by providing a uniqueness result when the coefficient is strictly increasing and 
C 1 . In that case, we show that the reconstruction amounts to of the inversion 
of the Laplace Transform of the Hankel Transform of a target function. We have 
not attempted to address the numerical resolution of such a notoriously ill-posed 
problem: it is very possible that it is manageable when restricted to a particular 
class of functions. 
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Appendix A. Proof of Proposition [2] 

Proposition 18. 14 For any two solutions Z\ and Z 2 of ([T]). Z T DZ 2 is inde- 

10 



pendent of t, where D = 

Proof. We compute that 
d 



10 
1 
10 



dt (zlDZ 2 ) =i(Zl)DZ 2 + ZlDf t (Z 2 ) 

= -\Z T B T DZ 2 + Z T D\BZ 2 
= i^Z T {DB-B T D)Z 2 , 
and since DB = B T D, we obtain that Z T DZ 2 does not depend on t. □ 

We now prove Proposition [2j 
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Proof of Proposition [H Since the entries of the matrix B is are piecewise continu- 
ous functions, global existence and uniqueness of the fundamental solution of the 
system is well-known see e.g. [15] . Now we prove the existence and uniqueness of 
solutions to the transmission problem via the fundamental solution. Let us define 
the fundamental solution of equation ([7]) as follows 

$(*) = (01, 02, 03, 04), 



where 0,- = 



nj,<p2j,<P3j,(P4j 



f , j = 1,2, 3, 4 with $(0) =Id. 



We have $(1)Z(0) = Z(l) therefore 



(41) 



where 



Ai(l) 



T 2 



M{1) 



o 



h 

h 



Ax 



and 





012 " 


- 011 


014 - 013 


-1 









022 - 


- 021 


024 - 023 


-1 









032 - 


- 031 


034 - 033 





-1 






042 " 


' 041 


044 - 043 





-1 






L 


-(011 + 012) 


-(013 


+ 014) 


1 





"(021 + 022) 


-(023 


+ 024) 





-1 


~(031 +032) 


-(033 


+ 034) 





1 


-(041 + 042) 


"(043 


+ 044 


) 



Ao 



To establish the announced existence and uniqueness property, let us now show 



that for any t € [0, 1], |detAi| > 2. Let U 



U 2 

v A 



andF 



Vi 

v 2 
v 3 

Va 



be given by 



Set Ci 





U x 


-u 2 




v 3 


-Vi 


" Vi- 


v 2 ' 


v 3 - 


Vi 



v x -v 2 
v 3 ~v A 



. We compute that 



Vi = i2 - 0a and Vi = i4 - 0i3- 
Then the determinant of A\ satisfies 

det(A)= % 
Ui-U 2 ] An \Vi- 

d-C 2 = (UiVi + U 2 V 2 - U X V 2 - U 2 Vi) + {U 3 V 3 + UiVi - U 3 V 4 - UiV 3 ). 

From the identity $(0) = Id , we deduce that V T {0)DU(0) = 0, Since U and V 
are two solutions of ([7]), Proposition [T51 shows that this implies V T DU — for all 
t € [0,1]. Expressing this identity in terms of the components of U and V, we 
obtain 

(XJ X V 2 + XJ 2 V X ) + (C/3V4 + C/4V3) = 0, 
and in turn 



(42) C x -C 2 = UiVi + U 2 V 2 + U 3 V 3 + UiVi 

Cauchy-Schwarz inequality thus shows that 

C x -C 2 = XN{U)N(V), 



U-V. 
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with A € C and |A| < 1. Similarly, starting from the identity U T _(0)DU(0) = -2 
and V T (0)DV(0) = -2 we obtain C x ■ Ci = N{U) 2 + 2 and C 2 • C 2 = JV(U) 2 + 2. 
Therefore 

Ci • Ci 02 
(7i • 02 C2 
|A| 2 )iV(£/) 2 iv (U) 2 + 4 + 2JV([/) 2 - 
>4. 

D 



|det(i4i)| s 



det(Ai )det(Ai) 



(1 



C x 
C 2 

>r2N{Vf 



Problem 
condition 

find r £ 



Appendix B. Proof of Proposition \W\ and Proposition [J] 

can be formulated slightly more generally as follows. Given an initial 



u(0) = l + r,v(0) = 1-r, 
t € C such that u(l) = v(l) = t, when u and v satisfy 

-u = if(x)v, —v = ig(x)u, 



dx ' dx 

with / and g two positive functions, bounded above and below by positive constants. 
This applies to system (JT5|) with v = Z3, u = z&, r/f = A, and rjXg = A 2 + £22 — 1- 
Note u is a solution of the real valued elliptic equations 

d (1 d \ 

(43) Tx\lTx w ) +9W = ^ 

The general uniqueness results for ordinary differential equations show that u is 
a linear combination of the two fundamental solutions wi and W2, of (1431) cor- 



f 



responding to the initial conditions wi(0) 
^w 2 (0) = 1. 

u(t) = Awi(t) +Bw 2 (t) 
Note that the Wronskian of the problem is constant: 



-4 

U ' /(0) dx 



/(0) dx 



wi(0) = and w 2 (0) 



/(■> 



wi— w 2 
ax 



W 2 — Wi 

dx 



1. 



Using the initial and final conditions, we obtain 



.4 
B 



1- 
ifl 



r 
-r) 



and 



wi(l) 



f(i)dx Wl W /(1 .'-;.' 
which leads to the following formulae for t and r 



w 2 (l) 
' &w a (l) 



.4 
B 



t 
if 



t 



e 2 (l)+iei(l) 



and r = 



iei(l)-e 2 (l) 
e 2 (l) + iei(l)' 



where e; 



d_ 

f(x) dx 



W{ — iwj for i = 1,2. Using the Wronskian identity, we note that 



t is always well defined, as a simple calculation shows that 
4 
W 



|ei(l)| 2 



e 2 (l)h 

2 



(tttT^W) +(/TT)^ Wl ( 1 )) +(w 2 (1)) 2 + (w 1 (1)) 2 + 2 



One can check that Itl 2 



= 1, therefore the amplitude of the reflected field, if 



it was available, would be redundant. 
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To prove that t is unchanged if / and g are replaced by f(x) = /(l — x) and 
g(x) = .9(1 — x), note the solution u corresponding to / and g can be written 
as linear combination of Wi(l — x) and 1^2(1 — x). The associated reflection and 
transmission coefficients f and t satisfy 



A 
-B 


= 


t 
it 


and 



w 2 (l) 



7TITSW2(1) 



/(i) 



A 

B 



1- 



wi(l) 

which shows that t = t and r = r. 
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